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|  ABSTRACT 

Hydrodvnamic  variables  such  as  density  ratio,  pressure  ratio, 
temperature,  fract'or.  of  dissociation  and  ionization  and  the  reflected  shock 
velocity  have  been  calculated  for  a  reflected  shock  wave  In  hydrogen  for  one 
dimension.  The  data  are  presented  as  families  of  curves  of  these  variables 
plotted  agarnst  the  initial  shock  Mach  number  for  a  range  of  unity  to  200  with 
the  initial  pressure  as  a  parameter 

Hydrodynamic  equations  are  derived  which  include  the  effects  o £ 
dissociation,  ionization  and  atomic  excitation  Reference  is  made  to  a  pre¬ 
vious  report  for  the  hvdrodynarmc  var'aV.e*  beh'nd  the  initial  shock  wave  in 
hydrogen  which  were  used  tn  these  calculations, 

Deviation i  from  the  so  ca'led  *  strong  shock  approximation" 
are  quite  large,  especially  <n  the  regions  where  the  gas  is  be4ng  dissociated 
or  ionized  The  density  and  pressure  ratios  across  the  reflected  shock 
wave  may  be  two  to  three  t»mes  greater  with  *  corresponding  decrease  in  the 
reflected  shock  vel<*c:*v  and  the  i em, 


erature. 
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1.  INTRODUCTION 

The  hydrodynamic  variables  for  a  primary  shock  wave  moving  into  hydrogen 

gas  at  room  temperature  (300*  K)  have  been  calculated  previously  and  the  results 

were  published  in  a  Space  Technology  Laboratories  report*.  These  data  allow 

one  to  determine  quickly  such  quantities  as  the  degree  of  ionization,  the  density 

* 

ratio,  the  temperature,  etc.  ,  from  a  measurement  of  the  shock  velocity  . 

Very  often  one  wishes  also  to  know  the  ariables  behind  the  reflected  shoe': 
at  the  end  of  a  shock  tube.  Therefore,  usir.fi  the  result®  of  the  above -mentioned 
report,  the  calculations  have  been  extended  to  include  the  hydrodynamic  variables 
behind  the  reflected  shock.  The  variables  calculated  were  the  pressure  ratio 
across  the  reflected  shock,  the  density  ratio,  the  temperature,  the  fraction 
dissociated  or  ionized  and  the  r*  lb  t  red  shock  velocity.  The  latter  is  very  useful 
because  the  theoretical  reflv<  W*d  shock  velocity  car.  be  compared  with  that 
measured  experimentally  to  determine  whether  the  gas  is  behaving  according 
to  theory  •-  in  particular,  whether  equilibrium  has  been  established. 

The!  v  cult  illations  are  for  one  -dimensional  flow  only,  and.  as  imph/d  in 
the  title,  equilibrium  has  been  assumed  throughout.  The  range  covoreu  in  this 
report  is  from. a  shock  Mach  number,  of  unity  up  to  200.  Above  M  -  2*>0 

th.  wi/iabivs  are  given  wnh  sufficient  accuracy  by  asymptotic  values. 


Turner,  hi.  B.  ,  “Equilibrium  Hydrodynamic  Variables  Behind  a  Normal 
Shod:  Wave  in  Hydrogen,"  STL  report  GM-TR -0 l65-004t>0,  dated 
lb  August  1953 

rt 

The  hydrodynamic  curves  given  m  GM-TR -0  165 -00460  and  in  the  present 
report  apply  also  to  the  isotopes  of  hydrogen,  namely,  deuterium  and  trituim. 
The  shock  velocity  must  therefore  be  measured  by  a  dimensionless  variable, 
M ,  which  is  the  shock  velocity  divided  by  the  sound  speed  of  the  undisturbed 

ga::. 
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II.  THEORY 

The  calculations  naturally  divide  into  several  regions  depending  upon  the 
conditions  of  the  gases  behind  the  primary  and  reflected  shock  wave*.  The 
formulae  and  methods  of  calculation  are  different  for  the  various  regions.  The 
regions  are  the  following: 

(a)  The  gas  behind  the  reflected  shock  can  be  considered  ideal,  Mg  =  1-5. 

(b)  The  gas  behind  the  reflected  shock  is  partially  dissociated,  Mg  =  4  -  12. 

(c)  The  gas  behind  the  reflected  shock  is  completely  dissociated  and  the 

gas  behind  the  primary  shock  is  partially  dissociated,  Mg  =  12  -  18. 

(d)  The  gas  behind  the  reflected  shock  is  partially  ionized  and  the  gas 
behind  the  primary  shock  is  not  completely  dissociated,  Mg  =  16  -  18. 

(e)  The  gas  behind  the  reflected  shock  is  partially  ionized,  Mg  =  20  -  44. 

(f)  The  gas  behind  the  reflected  shock  is  fully  ionized  and  the  gas  behind 

the  primary  shock  is  still  partially  ionized.  =  40  -  64. 

(g)  The  gases  behind  both  the  reflected  and  primary  shock  waves  are 
completely  ionized,  Mg  =  52  -  200. 

(h)  Above  Mg  =  200  asymptotic  values  can  be  used  with  sufficient 
accuracy. 

There  will,  of  course,  be  some  overlapping  of  these  regions.  The  notation 
used  is  as  follows : 

p  =  pressure 
p  =  density 
T  =  temperature 

a  =  fraction  dissociated  or  ionized,  as  the  case  may  be 
u  =  flow  velocity 
a  -  sound  speed 
U  =  primary  shock  velocity 
=  reflected  shock  velocity 

h  =  specific  enthalpy 
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Subscripts  are  used  to  designate  a  particular  region  of  the  flow: 

(1)  -  the  undisturbed  gas  at  room  temperature 

(2)  -  the  gas  behind  the  primary  shock  wave 

(5)  -  the  gas  behind  the  reflected  shock  wave 

Consider  now  a  primary  shock  wave  with  velocity  U  hitting  a  solid  wall 
and  reflecting  with  a  velocity  . 


Fig  1 

The  condition  for  rejection  is  that  the  gas  behind  the  primary  shock  wave, 
moving  with  a  velocity  u^  ,  be  brought  to  rest,  i.e.,  u^  =  0. 

One  method  of  approach  is  to  consider  the  more  general  problem  of  a 
shock  wave  moving  in  a  tube  with  flow  velocities  on  either  side  of  the  shock 
wave.  Let  the  region  ahead  of  the  shock  be  denoted  by  the  subscript  (a)  and 
the  region  behind  the  shock  by  (b). 
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Then  the  three  conservation  equations  across  the  shock  wave  can  be  written: 

(U  -  u^)  =  pa  (U  -  u^)  --  conservation  o£  mass,  (1) 

(U  -  -  pa  (U  -  ua)^  =  Pa  -  -  *  conservation  of  momentum,  (2) 

1  (U  -ua)2  +  ha  =  j  (U  -  --  conservation  of  energy .  (3) 

From  these  three  equations  can  be  derived  the  following  relation  for  the  flow 
velocities , 


Now  let  region  (b)  be  region  (5)  and  (a)  be  (2).  Then 


To  proceed  further  it  is  necessary  to  spec  if y  the  condition  of  the  gases 

on  either  side  of  the  reflected  shock  wave.  Consider  first  an  ideal  gas.  Then 

/  + 
the  ratio  P^/P^  18  g*ven  by  the  Rankine -Hugoniot  equation  . 


See  Reference  1. 
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where  p  is  given  by  =  y  +  l/y  -  1  .  The  ratio  p^/p^  can  written  as 


using  the  Rankine -Hugoniot  equation  for  p  /p.  and  the  expression  for  the 
Z  fc  1 

sound  speed,  a  =  y  (p/p)-  The  flow  velocity  behind  the  primary  shock  wave 
in  an  ideal  gas  is 


Z 


This  equation  follows  from  equation  (4)  and  equation  (6)  applied  to  the  primary 
shock  wave. 

Putting  these  into  equation  (3)  gives  the  quadratic  equation, 


which  has  only  one  physically  possible  root, 


(9) 
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The  density  ratio  across  the  reflected  shock  is  given  by  equation  (6),  and  the 
temperature  ratio  is  merely  the  pressure  ratio,  p^/p^,  divided  by  the  density 
ratio,  p^/p ^  From  the  conservation  of  mass  across  the  reflected  shock  it 
follows  that 


(10) 


There  are  two  ranges  of  shock  velocity  where  the  gas  may  be  considered 
ideal.  The  first  is  the  range  of  relatively  weak  shocks  before  the  hydrogen 
starts  dissociating.  Then 


v 


6  . 


The  second  is  the  range  where  the  gases  behind  both  the  reflected  shock  wave 
and  the  primary  shock  wave  are  completely  ionized  and  the  enthalpy  of  gas 
behind  the  primary  shock  is  much  larger  than  the  energy  of  ionization.  This 
range  of  shock  velocities  is  from  =  200  to  o o.  For  this  case 
y  =  5/3  so  |jl  =  4  .  Then 


and 


From  equation  (10) 


U2  3 

TT  ~  T  ’ 


1 

7 


but  since 
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The  temperature 


These  are  the  asymptotic  values  which  are  approached  for  very  high  shock 

velocities,  and,  as  wilt  be  seen  from  the  curves,  they  are,  in  most  cases, 

sufficiently  accurate  above  M  =  200. 

s 

Next  the  calculations  will  be  discussed  for  the  regions  where  the  gas  behind 
the  reflected  shock  wave  is  fully  ionized  or  fully  dissociated.  The  gas  behind 
t.he  primary  shock  can  be  partially  or  fully  ionized  or  dissociated,  respectively. 
This  mrludes  regions  (c),  (f)  and  (g)  as  outlined  previously. 

At  this  point  it  is  necessary  to  consider  the  equation  of  state.  It  can  be 
written  in  general  as 

p  =  XpRT  (11) 

For  an  ideal  gas  that  i.s  not  yet  dissociated  X  =  1.  For  a  partially  dissociated 
gas  X  =  (1  t  c)  where  a  is  the  fraction  dissociated,  and  for  a  partially  ionized 
gas  X  -2(1+  q  )  where  a  is  now  the  fraction  ionized. 

We  will  again  make  use  of  equation  (5).  The  density  ratio  across  the 
reflected  shock  is  now  given  by  a  revised  Rankine -Hugoniot  equation, 


4 


+  4-2 


1  -°2  6 

l  +  a7i  T2 


(U) 


g;  is  the  fraction  ionized  or  dissociated  and  0  is  the  energy  per  atom  of 
ionization  or  the  energy  per  molecule  of  dissociation  respectively  divided  by 
the  Boltzmann  constant. 


« 

Strictly  one  should  subtract  from  k0  the  average  excitation  energy  per  atom 
or  per  molecule  in  state  ( Z ).  This  was  done  for  the  case  of  partial  ionization, 
but  not  for  the  case  of  partial  dissociation  where  molecular  excitation  data 
was  not  available. 
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In  order  to  find  p^/p^  we  note  that 


i  pJl  ?J.  -  L  ?-L  J. 

X  p2  P1  ‘  X  p2  5a 


*Z  5\  *  Z  Z 

*1  =  “  TT  ‘  ' 


Putting  these  into  equation  (5)  results  in  the  quadratic  equation, 


2 

P-  P- 

—  -  2I2A+1-B)  —  -  (A+2B-1)  =  0, 

P2  p2 


(13) 


whe  re 


and 


A  -  7  ^  Tl 

*  -  nr  -j  tj 


B  -  I  _2- 

3  1  *a  i!  Ta 


The  values  of  variables  such  as  T^/Tj,  u^/a^,  and  were  taken  from 

Reference  1  After  solving  this  quadratic  equation  for  P^/p^»  p-/p^  can  be 

found  from  equation  (12),  The  reflected  shock  velocity  is  given  by  equation  (10). 

The  calculation  of  the  variables  behind  the  reflected  shock  for  shock 
velocity  ranges  (b) ,  (d)  and  (e)  was  somewhat  more  difficult  as  it  involved  an 
interpolation  procedure.  These  are  the  ranges  for  which  the  gas  behind  the 
reflected  shock  is  either  partially  ionized  or  partially  dissociated.  It  is 
necessary  to  assume  the  state  of  the  ga^  behind  the  reflected  shock  and  then 
calculate  the  residual  velocity  Two  or  three  calculations  are  usually  sufficient 
to  give  accurate  interpolated  values  for  the  variables  for  zero  velocity. 
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The  equation  for  the  \  elocity  for  the  case  of  partial  ionization  or 

* 

dissociation  was  derived  in  Reference  1  for  the  primary  shock  wave  ,  but 
with  a  slight  generalization  it  can  also  be  used  for  the  reflected  shock  wave. 
For  the  reflected  shod;  wave  the  equation  can  be  written 


uc  -  u  \ 
'5  Z  ; 


5x 

7  XZ  Tf 


KS  T. 

q  t~2 


S  T5 

*2  T2 


4<W 


±lRTU 


{hs-hi]  st5  ; 

\>RT2  *71  ‘  . 


The  procedure  is  to  first  choose  the  enthalpy,  h-  .  Then  the  pressure 
ratio  across  the  reflected  shoe  :  wave  is  determined  by  the  relation, 


which  is  derived  from  the  conservation  conditions  across  a  shock  wave.  From 
h^  and  p^  one  can  now  determine  and  .  This  can  be  facilitated  by 

using  a  family  of  curves  of  enthalpy  versus  temperature  for  several  pressures 
and  a  family  of  curves  of  fraction  ionized  or  dissociated  versus  temperature 
for  the  same  pressures.  The  calculation  of  the  enthalpy  and  fraction  of 
ionization  is  explained  in  Reference  1,  and  the  curves  are  shown  in  Figures  3and4 
of  this  report. 

Similar  curves  for  the  enthalpy,  together  with  curves  for  the  fraction  of 
dissociation,  are  given  in  the  Bureau  of  Standards  Circular  364,  1  November  1953, 
entitled  "Tables  of  Thermal  Properties  of  Cases."  These  curves  cover  the 
region  of  dissociation  while  the  curves  in  Figures  3  and  4  cover  the  region  of 

* 

This  equation  was  originally  derived  by  E.  L.  Resler  in  his  Ph.  D.  thesis 
entitled  "High  Temperature  Gases  Produced  by  Strong  Shock  Waves." 

Cornell  University,  1951. 
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temperatures  where  the  hydrogen  is  partially  ionized.  The  enthalpy  for  partially 
ionized  hydrogen  given  in  Figure  3  includes  the  energy  of  excitation  of  the 
neutral  atoms.  The  Bureau  of  Standards'  curves  for  the  enthalpy  of  partially 
dissociated  hydrogen,  however,  do  not  include  the  effect  of  molecular 
excitation  of  the  undissociated  molecules. 

The  values  of  the  variables  in  region  (2)  such  as  u^/a^,  T^,  an<*  ^2 

were  taken  from  the  first  hydrodynamic  report*.  Two  or  three  estimates  of 
h^  must  be  made  to  give  values  of  u^  on  either  side  of  zero.  Then  the 
hydrodynamic  variables  h^,  p $/p£>  T^  and  can  be  interpolated  for  u^  =  0, 
The  density  ratio  is  given  by 

Ps  _  P5X2TZ 

p2  '  *71?$  \ 

and  the  reflected  shock  velocity  is 


U 


R 


u 


2 


The  hydrodynamic  variables  behind  the  reflected  shock  wave,  calculated 
according  to  the  methods  given  above,  have  been  plotted  on  graphs  shown  in 
Figures  5-9.  In  Figure  5  the  density  ratio  behind  the  primary  shock  wave 
has  been  plotted  along  with  the  density  ratio  behind  the  reflected  shock  wave  so 
that  the  total  density  ratio  can  be  quickly  determined  without  reference  to  the 
first  hydrodynamic  report*.  In  Figure  8  the  temperature  behind  both  the 
reflected  and  primary  shock  waves  are  shown. 

These  hydrodynamic  curves  show  the  large  effect  that,  dissociation  and 
ionization  have  on  the  variables  behind  the  reflected  shock  wave.  An  assump¬ 
tion  of  the  asymptotic  value  of  2.5  for  the  density  ratio,  for  example,  would 
be  too  low  by  a  factor  of  more  than  2  over  the  range  of  Ms  =  6  to  40.  Similar 
large  deviations  occur  in  the  pressure  ratio  across  the  reflected  shock  wave. 
The  asymptotic  value  is  6,  but  the  pressure  ratio  has  maxima  of  more  than  20. 
As  one  should  expect,  the  reflected  shock  velocity  becomes  much  less  in  the 
range  where  the  pressure  and  density  ratios  are  large.  The  minimum  reflected 
shock  velocity  is  less  than  one -third  of  the  asymptotic  value  of  l/l  U. 
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Specific  enthalpy  of  hydrogen  vs  temperature  for  several  pressur 
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10.000  20,000  30,000  40,000  50,000 

TEMPERATURE.  °K 

Fig-  4  Fraction  of  hydrogen  ionizeu  vs  temperature  for  several  pressure 
4c*1rulau*'J  from  the  Sriha  eutialinu) . 
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PRIMARY  SHOCK  VELOCITY. 
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Fig.  7  Fraction  ionized  or  dissociated  behind  the*  reflected  shock  wave. 
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PRIMARY  SHOCK  VELOCITY. 


